Effect of uniaxial strain on plasmon excitations in 
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Abstract. Uniaxial strain is known to modify significantly the electronic properties of 
graphene, a carbon single layer of atomic width. Here, we study the effect of applied strain on 
the composite excitations arising from the coupling of charge carriers and plasmons in graphene, 
i.e. the plasmarons. Specifically, we predict that the plasmaron energy dispersion, which has 
been recently observed experimentally in unstrained graphene, is shifted and broadened by 
applied uniaxial strain. Thus, strain constitutes an additional parameter which may be useful 
to tune graphene properties in plasmaronic devices. 



Graphene is a monoatomic carbon layer with honeycomb structure. After having been 
considered for decades as the ideal constituent of most compounds of carbon in the sp 2 
hybridization state, it has been recently obtained in the laboratory [lj, thus kindling intense 
research activity both on the experimental and on the theoretical side. Graphene is especially 
characterized by a quasiparticle band structure consisting of two bands, touching at the Fermi 
level in a linear, cone-like fashion at the so-called Dirac points ±K, and a linearly vanishing 
density of states (DOS) at the Fermi level [2j [3]. These peculiar electronic properties, along 
with the reduced dimensionality, have remarkable effects on the electromagnetic properties 
of graphene. These include, e.g., the reflectivity [4], the optical conductivity [5j El El [9], 
the plasmon dispersion relation [10], [Til EZ1 03], as well as a newly predicted transverse 
electromagnetic mode [2], which is characteristic of a 2D system with a double band structure, 
such as graphene. 

The composite elementary excitations arising from the coupling of charge carriers and 
plasmons, the so-called plasmarons, have been considered in a general context earlier on by 
Lundqvist |15l I16| . Recently, plasmarons have been experimentally observed in graphene by 
means of angular resolved photoemission spectroscopy (ARPES) [T7|, and their dispersion 
relation described theoretically within an improved version of the random phase approximation 
(GoW-RPA) [TB]. 

In n doped graphene, a plasmaron mode with momentum k results from the relatively strong 
coupling of a quasihole with momentum k+q, and a plasmon with momentum — q, the quasihole- 
plasmon coupling being stronger when the two excitations have the same group velocity |18j . 



At k = 0, the plasmaron relative momentum modulus turns out to be 
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where [i is the chemical potential, vf the Fermi velocity, e = eoe r the dielectric constant. 
Therefore, the plasmaron binding energy with respect to the Fermi energy can be estimated, in 
first approximation, as the sum of the energies of the bare quasihole and plasmon, both having 
momentum modulus q, viz. 

E P = -»-a-^, (2) 
vf 2e r 

where a is the fine structure constant. In the realistic case of graphene on a SiC>2 substrate, 
Eq. ([5]) yields Ep ~ — 1.25/i. A more accurate estimate, including the contribution of the 
quasihole-plasmon interaction at the GoW-RPA level [18], yields Ep ~ — 1.3/x, in better 
agreement with the experimental results [17]. 

We now consider the effect of uniaxial strain on the graphene sheet. It has been shown that 
this amounts to a shift of the Dirac points in reciprocal space, and to an anisotropic deformation 
of the Dirac cone centred at those points |19j . with elliptic sections at constant energy. Such 
a strain-induced angular dependence may also be nonuniform, as in the case of coordinate- 
dependent strain [20] . The effect of applied uniaxial strain on the plasmon dispersion relation 
of graphene has been studied in Refs. [12], [13]. Besides, is has been shown [19] that strain 
may significantly modify the dispersion relation of a transverse plasmon mode, which has been 
recently predicted to occur in graphene [T3]. Therefore, it can be expected, on general grounds, 
that strain affects the energy dispersion of the plasmaronic modes. 

Denoting by e the modulus of applied strain, 9 its direction (with 8 = and 6 = ir/6 referring 
to strain along the zig zag and armchair directions, respectively), and v the Poisson's ratio of 
graphene {y = 0.14 [21], to be compared with the known experimental value v = 0.165 for 
graphite [22], and with v = — 1, corresponding to the hydrostatic limit), following Ref. [19], the 
quasiparticle dispersion relation, to linear order in e, reads 



e q = ±hv F q[{l - k(1 - v)e) - k(1 + ^)ecos(26» + 2cp)], (3) 

where the + (— ) sign refers to the conduction (valence) band, and while the plasmon dispersion 
relation under strain becomes 



frup\ = y [1 - k(1 + v)e cos(20 - 2y>)] y/q. (4) 

Deriving the corresponding group velocities from the above Eqs. (|3|) and (HJ), Eqs. (pQ) and ([2]) 
for the plasmaron momentum and energy, respectively, get modified into 
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to linear order in the strain modulus e. Here, k = (a/2t)\dt/da\ — \ ~ 1.1 is related to 
the logarithmic derivative of the nearest-neighbor hopping t at e = 0, a is the carbon-carbon 
distance, and 9 is the direction of the stress. 

Eq. (|5b|) shows that, in the presence of applied uniaxial strain, the plasmaronic energy at 
k = acquires an explicit dependence on the angle ip of the quasihole momentum q. This is 
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due to the anisotropy of both the electronic and the plasmon spectrum. Correspondingly, the 
plasmaron energy is characterized by a central value 

E c P = -fi-a—^-[l + K(l-iy)e], (6) 
vf 2e r 

and a strain-induced energy spread 

AE P = a— —n(l + u)e. (7) 

Considering again the realistic case of graphene on a Si02 substrate, one can estimate the 
central plasmaron energy in the unstrained case as Ep(e = 0) = —1.25// = —125 meV, for 
[i = lOOmeV, with zero energy spread. Correspondingly, in the case of an applied strain 
e = 10 %, one finds a central plasmaron energy of Ep(e = 10 %) = —127.4 meV, with an 
energy spread AEp(e = 10 %) = 6.27 meV. 

In conclusion, the effect of applied uniaxial strain on graphene is therefore that of shifting and 
broadening the plasmaron energy, proportionally to the strain modulus. Therefore, by suitably 
applying uniaxial strain, one gains further control on the energy of the plasmaronic excitation, 
besides the possibility of tuning the relative dielectric constant e r |23| . This may be instrumental 
for the realization of 'plasmaronic' devices. 
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